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I ■ ABSTRACT 

O I We use numerical simulations to investigate how the statistical properties of dark 

matter (DM) haloes are affected by the baryonic processes associated with galaxy for- 
mation. We focus on how these processes influence the spin and shape of a large number 
of DM haloes covering a wide range of mass scales, from dwarf galaxies to clusters, 
at redshifts zero, one and two. The haloes are extracted from the Overwhelmingly 
Large Simulations, a suite of state-of-the-art high-resolution cosmological simulations 
^ ■ run with a range of feedback prescriptions. We find that the median spin parameter 

in DM-only simulations is independent of mass, redshift and cosmology. Baryons in- 
crease the spin of the DM in the central region (< 0.25 r2oo) by up to 50 per cent when 
feedback is weak or absent. This increase can be attributed to the transfer of angular 
momentum from baryons to the DM. We also present fits to the mass dependence of 
, the DM halo shape at both low and high redshift. At ^ = the sphericity (triaxiality) 

is negatively (positively) correlated with halo mass and both results are independent of 
cosmology. Interestingly, these mass-dependent trends are markedly weaker at z = 2. 
While the cooling of baryons acts to make the overall DM halo more spherical, stronger 
feedback prescriptions (e.g. from active galactic nuclei) tend to reduce the impact of 
baryons by reducing the central halo mass concentration. More generally, we demon- 
strate a strongly positive (negative) correlation between halo sphericity (triaxiality) 
■ and galaxy formation efficiency, with the latter measured using the central halo baryon 

fraction. In conclusion, our results suggest that the effects of baryons on the DM halo 
spin and shape are minor when the effects of cooling are mitigated, as required by 
realistic models of galaxy formation, although they remain significant for the inner 
halo. 

Key words: galaxies: haloes - galaxies: clusters: general - galaxies: evolution - meth- 
ods: numerical - cosmology: theory 
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1 INTRODUCTION 

A natural consequence of the standard hierarchical struc- 
ture formation paradigm is that the shapes of dark matter 
haloes are triaxial, a property t hat is inherited from their 
progenitor density perturbations (|Bardeen et al.|[l986l V This 
additionally leads to aspherical growth as the halo accretes 
matter from preferential directions, associated with the sur- 

* E-mail:sarah. bryan@gmail.com 



rounding sheets and filaments. The anisotropic accretion his- 
tory of a halo also affects its angular momentum distribu- 
tion, through the presence of non-zero torques. It is therefore 
clear that both the shape and spin of a dark matter halo are 
important diagnostics for an accurate determination of their 
structure and formation history. 

While the spin and shape of a dark matter halo 
are not directly observable, they have important conse- 
quences for the structure and dynamics of galaxies. For 
example, halo spin is an important parameter in galaxy 
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formation models as it affects the size of t he embed- 
ded galactic disc (e.g. I Mo et all 1 19981 . but see I Sales et al] 
I2OI2I ). Deviations from axisymmetry in elliptical galax- 
ies a re likely to influence th e gas kinematics of the sys- 
tem (|de Zeeuw fc Franxlll989l ). and may be responsible for 
exciting or s ustaining warps and stabilisin g or deforming 
polar rings (jSteiman-Cameron et al.l Il992l ). Axisymmetry 
may also influe nce the fuell ing efflciency of the central 
black hole ((Franx et al.l Il99ll '). Misalignment of the angu- 
lar momentum of the halo and the galaxy may be re- 
sponsible for the anisotropic distri butio n of subhaloes and 



satellite galaxies (i Holmbera ll96E and iKnebe et al 



and could ca use galac tic warps I Ostriker fc Binnev 



2004) 



iDebattista fc Sellwoodill999l : lBailin fc Steinmet j|2004h 



198£ 



On cluster scales, asphericity in the dark matter halo 
will naturally correspond to asphericity in the gas den- 
sity and will impact the shape of X-ray isophotes and 
the Sunyaev-Zel'dovich signal. Understanding the intrin- 
sic shape of dark matter haloes is also important for 
weak lensing analysis (see, for exam ple, the discussions in 
iBecker fc Kravtsovll201ll : iBettI 120121 ') and it is well known 
that intrinsic ellipticity can con tribute signiflcantl y to a lens- 
ing halo's ability to form arcs (|Oguri et al.ll2003l ). 

Several methods are used to constrai n galaxy and 
halo shapes obse rvationally (see, for example. ISackettll 19991 : 
iMerrifleldl 120041 '). Unfortunately, different observational 
methods currently yield systemat ically different results for 
the shapes of dark m atter haloes (|011ing fc Merrifleldll2000l : 
lO'Brien et al]l2010l '). However, it is worth noting that the 
observations cover a large range of systems and vary in the 
extent to which the halo is probed, making a direct compar- 
ison somewhat difficult. Whether the discrepancies result 
from halo-to-halo scatter or from systematic errors in the 
observed estimates is unclear. Given the rapidly accumulat- 
ing data set and ever increasing sophistication of the data 
analysis tools, one can soon expect to have an observational 
data set which may be more directly compared with theory. 

Theoretically, the predictions for the distribution of 
angular momentum and halo shapes have been studied ex- 
tensively, primarily using cosmological iV-body simulations 
dFrenk et al]|l 988': 'Dubinski fc Carlborg"l99l'; ' Warreii et all 
19921: ICole fc Lacey.,199& : .Bullock. 2002: Jing fc Sutdl2002l: 



Bailin fc Steinmet j 200E : lAUgood et al.l20Q6l : iBett et all 



20071: iMaccio et al' '2008*: Jccson-Danicl ct al.' '201 1': 



'Vera-Ciro et al. 201 1: Bctt 2012; Zcmp ct al. 2012). There 



is a general consensus that cold dark matter haloes have 
approximately log-normal spin distributions and are triax- 
ial, with sphericities, (c/a) ~ 0.5 — 0.8 Q and elongations, 
(6/a) ~ 0.4 — 1. Furthermore, haloes are generally found 
to be highly flattened and show a tendency toward prolate 
shapes [c/h > h/a), especially in the inner regions. There is 
also general agreement that the sphericity decreases with 
increa sing halo mass and that the spin is independent of 
mass (lBuUock ''2002': ji iig fc Suto '2002: Springel et al. '2004*; 



Hopkins et al.. 2005; .Bett et al.l ,2007 : Maccio et al., ,200a: 



Jeeson-Daniel et al.ll201lh . In addition. Ijeeson-Daniel et alj 



l|201ll ) have shown that sphericity is strongly correlated 
with concentration, while both triaxiality and spin are 
anti-correlated with concentration. 



^ Where a > fe > c are the eigenvalues of the halo's inertia tensor. 



While these results are interesting, a significant 
uncertainty is how the dark matter halo is affected by 
the additional, non-gravitational processes acting on the 
baryons. Recent work has clearly established that the 
condensation of baryons to the centre of dark matter 
haloes tends to increase the cent ral angular momen- 



tum of the halo (see for e xample, Sharma fc Steinmetz 
20051 : iKaufmann et all l2007l : lAbadi et al.ll2010l : iBett et al 



2010|) and to make the h alo more spherical or axisym- 



metr ic (see, for example. iKatz fc Gumil 199ll: Dubinski! 



19941: lE^ ard et al.' ' 19941; iBarnes fc Hernguistl Il996l: 



Tissera fc D omingucz-Tcnreiri 
2004 



Pedrosa et al 



19981: iSprineel et al l 



Kazantzidis et al.l 120041: Dcbatti sta et aL. 
2OI0I : iTissera et al.li201Q : .Brvan et al.l 



2008; 



20121 : 



Zemp et al.ll2012l '). This result has been used to explain 



the discrepancy between the strongly prolate-triaxial shape 
found in A'^-body simulations and the more spherical 
observed systems. 

Incorporating baryonic physics in cosmological simula- 
tions is a non-trivial task and the computational cost of this 
process has placed limits on both the parameter space and 
the size of the sample of haloes explored to date. The de- 
tailed nature of the baryonic processes involved in galaxy 
formation and the precise influence of these processes on 
galaxies therefore remains largely uncertain. In this paper, 
we attempt to make progress on both fronts, by studying the 
spin and shape distributions for a large (> 1000) sample of 
dark matter haloes, spanning a range of mass (from dwarf 
galaxies to clusters) and redshift (z = 0, 1 and 2). We do 
this using the Ove rwhelmingly Large Simulations (OWLS; 
IS chave et al.ll2O10l ) - a suite of cosmological hydrodynamical 
simulations run with many different physical prescriptions 
for the baryons. By providing identical simulations run with 
different implementations of the subgrid physics, OWLS of- 
fers the opportunity to explore the effects of baryons under 
a range of physical conditions, for the same population of 
haloes. In particular, we use the OWLS data to quantify, in 
a statistically meaningful way, the influence of feedback pro- 
cesses (from no feedback, to feedback from stars and black 
holes) on the spin and shape distributions of dark matter 
haloes. 

The paper is organised as follows. The simulations used 
in this analysis, and the methods used to identify haloes and 
to estimate their spins and shapes are outlined in section (2] 
Our results are presented in section 13.11 including fltting 
formulae for the predicted correlations between halo shape 
and mass. The robustness of our results is demonstrated 
via a resolution study, given in the appendix. Finally, we 
summarise our main results in section 2] 



2 METHODOLOGY 
2.1 Simulation details 

The haloes used for this analysis were extracted from a sub- 
set of the OWLS runs. For detailed inf ormation abo u t thes e 
simulations the reader is referred to ISchave et al.l (|201(j ): 
here, the most relevant aspects are briefly reviewed for con- 
venience. 

For all runs, cosmological initial conditions were set 
up using a transfer function generated with CMBFAST 
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Table 1. A list of the OWLS runs used in this analysis. We use the same identifier for each run as in lSchave et al and comment 

on the subgrid physics implemented in each case. Key global properties of the simulations are also listed, namely the number of dark 
matter and baryonic particles; the comoving box length; and the dark matter particle mass. Values are presented for all runs analysed 
at ^ = 0. The 100/i~^Mpc boxes are also used for results at 2 = 1, while 25/i~^Mpc boxes are used at 2 = 2 (values for these runs are 
given in brackets). The maximum force softening was 0.5, 2.0 and 8.0 /i^^kpc for the 25, 100 and 400/i~^Mpc boxes respectively. 



Name 


Description 


A'"dm 


-^baryons 


Box length 


rriDM 










(/i-i Mpc) 


(fe-lMo) 


DMONLY 


Dark matter only runs 












WMAPl 


216^ 




50 


8.6 X 10** 




WMAP3 


5123 




100 (25) 


4.9 X 10** (7.7 X lO*') 




WMAP3 


5123 




400 


3.1 X 10"' 




WMAP5 


512^ 




100 (25) 


5.3 X 10** (8.3 X 10'^) 




WMAP5 


5123 




400 


3.4 X lOi" 


N0SN310ZC00L 


No feedback, primordial abundances for cooling 


5123 


5123 


100 (25) 


4.1 X 10** (6.3 X 10^) 


REF 


Weak stellar feedback, metal cooling 


5123 


5123 


100 (25) 


4.1 X 10** (6.3 X 10*^) 


WDENS 


Strong stellar feedback, metal cooling 


5123 


5123 


100 (25) 


4.1 X 10** (6.3 X 10^*) 


AGN 


Weak stellar & AGN feedback, metal cooling 


5123 


5123 


100 (25) 


4.1 X 10* (6.3 X 10^) 



l|Seliak &: Zaldarriaga|[l99^ ). Initial ( z = 127) positions and 
velocities were computed using the IZel'dovichl (|l970l 'l ap- 
proximation from an initial glass- like state ( Whitelll99 6l). All 
simulations we re run using a modified version of GADGET-3 
ijSpringellbOoBI ). For a given box size the same initial condi- 
tions were used in each run, allowing us to directly compare 
the same haloes evolved using different prescriptions for the 
sub-grid physics. A summary of the simulations used in this 
analysis, including values for key numerical parameters, can 
be found in Table [1] 

We begin by exploring dark matter only (DMONLY) 
simulations, as this allows us to validate our results by 
comparing to the existing literature, as well as to set 
the scene for exploring the impact of baryons. We con- 
sider three sets of dark matter only simulations, each run 
with different values for the cosmological parameters. Our 
main results (including baryons) assume values taken from 
the 3rd ye ar Wilkinson Microw ave Anisotropy Probe data 
( WMAPi: [Spergel et al.l l2007f) with [n™,nA, J^b, erg] = 
[0.238,0.762,0.0418,0.95,0.74]. This model was run with 
dark matter only, using 5123 particles in three box sizes. 
The two larger boxes (100 and 400 ft" ^ Mpc) were run to 
z — Q while a smaller, high-resolution box (25 /i~^Mpc) was 
run to 2: = 2; we use the latter to study the properties 
of haloes at high redshift. The maximum physical soften- 
ing length in the 25ft~^Mpc box was 0.5 /i~^kpc, while 
in the 100 and 400/i~^Mpc it is 4 and 16 times larger re- 
spectively. We also present resu lts for the newer WMAP5 
cosmology ijKomatsu et al.l|2009h with [0,mfiA,^b,n, as] = 
[0.258,0.742,0.0441,0.963,0.796] using the same boxes as 
before. Fin ally, we consider a run with the WMAPl 
cosmology l|Spergel et al.l [ioO^ ) with [Qrn,^A,^ib,n,as] = 
[0.25, 0.75, 0.045, 1, 0.9]. This simulation was run using 216^ 
particles in a 50 Mpc box, matching the resolution of 
the Millennium Simulation (Springcl et al. 2005), which was 
also run with the WMAPl parameters. 

To study the effect of varying levels of feedback on the 
spin and shape parameters of haloes extracted from ACDM 
simulations, four baryon runs from the OWLS simulations 
are considered, all run with the WMAP3 cosmology. These 
runs model the gaseous component using smoothed particle 
hydrodynamics (SPH). 



All of the sim ul ations include radiative 
cooling JWiersma et al | 2009al ). star formation 
(Schave & Dalla Vecchia] 200a ) and metal enrichment 
(Wicrsma et al. 2009b), but differ in their feed back pre- 
scriptions (for more details see iDalla Vecchia fc Schavd 
(200^). All models followed the timed release of H, He 
and 9 different heavy elements produced b y massive stars, 
AGB stars and supernovae of type I and II (|Wiersma et al.l 
|2009b) The first run (NOSN JSfOZCOOL) did not include 
any feedback processed; the second (REF) includes weak 
feedback from stellar winds and supernovae; the third 
(WDENS) includes stronger stellar feedback with a wind- 
speed that depends on the local gas density; and the final 
run (A GN), feedback from bo th stars and active galactic 
nuclei l|Booth fc Schavd \200^ ). In the feedback models, 
enhanced cooling from metal lines was also accounted for. 
All four models were run using the 100/i~^Mpc box (to 
z = 0) and the 25/i"^Mpc box (to z = 2). The same 
number of gas particles as dark matter particles (5123) 
adopted for each box. 



2.2 Halo sample 

Haloes were extracted from each sir nulation usi n g the 
Friends-of-Friends (FoF) algorithm (|Davis et al.l Il985f ) 
which links particles together within a fixed comoving sep- 
aration. This separation, known as the linking length, was 
set to 0.2 times the mean inter-particle distance. The FoF 
groups were then decompo sed into self-bound sub-haloes us- 
ing th e SUBFIND algorithm (|Springel et al ]|200ll ; lDolag" et al.l 
l2008h . Finally, a sphere was grown around the most bound 
particle of the most massive sub-halo until the mean total 
mass density was equal to 200 times the critical density. We 
choose to define our haloes as including all particles con- 
tained within this sphere. The mass and radius of the halo 
are referred to as M200 and r2oo respectively. 

Only haloes that contain at least 1000 particles are con- 
sidered in this analysis, as this ensures that the results are 



^ This simulation also neglects enhanced cooling through metal- 
lines. 
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fully converged (see the appendix for a discussion of the ef- 
fects of resolution). While estimates of the spin parameter 
are found to b e well resolved fo r haloes with more than 300 
particles (as in lBett et al.ll2007l ) a larger number of particles 
is required to resolve the halo shape (i n particular the triax - 
iality of the halo) , in agreement with iMaccio et al.l ()2008l l . 
This cut imposes a minimum halo mass of4.9xl0^^ /i~^M0 
at z = and 7.7 x 10^ h'^ Mq at z = 2, for our WMAP3 
runs. 

We have investigated whether our results are affected by 
restricting our sample to dynamically relaxed haloes. We es- 
timated the dynamical state using the centroid shift, defined 
to be the distance between the minimum potential position 
(the halo c entre) and the cent re of mass of the halo. In accor- 
dance with lNeto et al.l l|2007l ). a halo is defined to be relaxed 
if the centroid shift is less than 0.07 r2oo- We find that while 
relaxed haloes are typically more spherical and less triaxial, 
this restriction does not affect the trends found in this work 
and so we chose to present our analysis for the total halo 
sample (except where explicitly stated otherwise). Not only 
does this result in larger numbers of objects, our choice is 
also motivated by the fact that an observational cut based 
on the relaxation state of a halo is not a straightforward 
task. 

The numbers of haloes in our final samples are pre- 
sented in Table (2] Note that the number of haloes in 
DMONLY WMAP3 at z = is almost twice as large as 
in the baryon runs; this is due to the inclusion of haloes 
from the 400ft"^Mpc box in the former case. Overall, our 
DMONLY sample spans around three orders of magnitude 
in mass (from galaxy to cluster scales, i.e. M200 ~ lO^'^ — 
10^^ h~^MQ) at 2 = 0. Our z = 2 sample covers a similar 
dynamic range but at lower mass (dwarf galaxy to group 
scales, M200 ~ lO^'^ - 10^^ /i-^Mq). 

2.3 Defining iialo spin and shape 

The dimensionless spin parameter A provides a useful mea- 
sure of the amount of rotational support present within a 
dark matter halo. We estimate this property for the dark 
matter partic l es us ing the modified expression given by 
iBullock et all (|200ll ) 

A'(r) = , (1) 

where J is the angular momentum within a sphere of radius 
r, containing mass M and Vc = yX?M(<r)/r is the halo 
circular velocity at this radi us. This expre ssion reduces to 
the standard spin parameter l|Peebleslll96l ) when measured 
at the virial radius of a truncated singular isothermal halo. 
For a compari son of the dif f erent definitions for the spin 
parameter, see iMaccio et all l|2007l '): hereafter we drop the 
prime and refer to the modified spin parameter as A. 

To characterise the halo shape, we use the mass 
distribution tensor M. which has b een used exte n sively 
in the halo shape lite rature (e.g. ICole fc LacevI 1 19961 : 
iBailin fc Steinmet j|2005h . The components of the tensor (a 
square matrix) are 

Mij ^^^mkrk,trk,j, (2) 
fe 

where the index k runs over all dark matter particles within 



a given radius, mk is the mass of the fcth particle and rk,i the 
ith component of its position vector from the halo centre. 
The square roots of the eigenvalues of the mass distribution 
tensor, obtained using Jacobi transformations, are defined as 
a, b, c (where a > b > c) and are used to measure the shape 
of the simulated haloes. Note that the shapes obtained using 
the inertia tensor and the m ass distribution tensor are 
equivalent (|Bett et al.ll2007l ). 

Our shape results are presented using the following pa- 
rameters: s — c/ais used as a measure of halo sphericity; e = 
b/a as a measure of elongation; and T — {a? — b'^)/{a? — (?) 
as a measure of the triaxiality of the halo. A purely spherical 
halo will have s = e = 1 with T being undefined. Low values 
of T (T — >■ 0) correspond to oblate haloes while high values 
(T — > 1) correspond to prolate haloes. 

We note that computation of the mass tensor in a spher- 
ical region biases the shape towa rds higher sphericity; this 
is corrected for, as suggested in iBailin fc Steinmet j |2005| . 
by re-scaling the axis ratios s — >■ and e — >■ . While 
our adopted method described above may not be the most 
robust way of de scribing the phys ical halo sh ape (e . g. see 
the discussion in IZemp etliLlbOllI ). we follow iBettI (|2012l ) 
and use this simple approach as it is most directly compa- 
rable with observations and is adequate for the comparison 
we wish to present here. We checked that using the reduced 
mass distribution tensor, and also using an iterative tech- 
nique to determine the halo shapes, do not result in system- 
atically different results. 

Our main results for spin and shape parameters are pre- 
sented for the dark matter within r2oo, but we also consider 
results within the central region (0.25 r2oo) where baryonic 
physics plays a more signi ficant role. The convergence ra- 
dius fsee lPower et al.ll2003l ) is less than 0.25 7-200 for almost 
all of the haloes considered here. (The haloes that are not 
converged within 0.25 r2oo are excluded from the studies of 
the central region. In all cases this is less than one per cent 
of the sample) . Resolution tests for this region are discussed 
in the appendix. 



3 RESULTS 

We begin by presenting our results for the halo spin param- 
eter, then go on to explore the shapes of a large number of 
haloes selected from our cosmological sample. In both cases, 
we first discuss the results from dark matter only simula- 
tions before studying the impact of the baryons. Quantities 
are measured within the central region (0.25 r2oo) and at 

?'200 

3.1 Halo Spin 

3.1.1 Dark Matter Simulations 

The distribution of dark matter halo spin is typically found 
to be well characterised by a log- normal distribution 

^ Although we note that a more accura te fitting fo r mula for the 
distribution of halo spins is presented in lBett et al. I 1I2OO7I I. 
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Table 2. Number of haloes (containing at least 1000 particles) and median halo mass M200 (in lO^^/i ^Mq) for each simulation at 
^ = 0, 1 and 2. 



Simulation 


2 = 
-'^haloes 




M200 


z 
■ 


1 

M200 


z = 
-'^haloes 


2 

M200 


DMONLY WMAPl 


431 


1.7 










DMONLY WMAPZ 


7094 


2.3 


3765 


0.77 


5334 


0.030 


DMONLY WMAPb 


8188 


3.6 










AGN 


3878 


0.88 


3375 


0.83 


4743 


0.033 


WDENS 


3884 


0.91 


3371 


0.84 


4699 


0.035 


REF 


3990 


0.99 


3580 


0.93 


4832 


0.044 


NOSN.NOZCOOL 


4634 


1.0 


4041 


0.91 


5683 


0.040 



0.8 I 

0.7 - z = 
''200 



0.6 



Med M, 



0.5 - n = -3.34 
a = 0.62 



D- 0.4 - 
0.3 
0.2 - 
0.1 



J 



/ 



In [V) 



0.7 - z = 
0.25r. 

0.6 - 



200 

Med M2oo = 2x10^^/7'^ 



0.5 - n = -3.34 
o = 0.62 



Q- 0.4 - 
0.3 - 
0.2 
0.1 - 



Ot 



A 



/ 



-2 -1 



0.8 
0.7 
0.6 
0.5 

Q- 0.4 
0.3 
0.2 
0.1 




z = 2 

''200 



|i = -3.26 
= 0.6 



1 



J 



/ 



\ 



In (X) 




In (X) 



In {%) 

Figure 1. The distribution of log halo spin parameter for dark matter particles in haloes at 2 = (left) and 2 = 2 (right), taken from 
the DMONLYWMj4P3 runs. The top (bottom) panels show the spin distribution computed using all particles within r2oo (0.25r2oo)- 
The mean and standard deviation of the distribution are computed for each sample and are listed in Table |3] The best-fitting Gaussian 
curve, assuming these parameters, is overlaid. For comparison, we over-plot arrows representing the mean of the distributions for the two 
baryon runs with the most extreme feedback prescriptions - the run with no feedback (NOSN_NOZCOOL; green arrow) and the AGN 
feedback run (AGN; red arrow). It is evident from this figure that efficient cooling results in an increased DM spin parameter, especially 
within the central region of the halo, while strong (AGN) feedback results in values that are indistinguishable from the dark matter only 
case. 



where Ao and a are free parameters, determining the 
mean and standard deviation respectively. We confirm 
this result in Fig. [Jl where the distribution from the 
DMONLY WM^P 3 haloes is shown as a histogram. Values 
for the mean and standard deviation of the sample are also 
given, and are used together with Equation [S] to predict the 
equivalent log-normal spin distribution (shown as the solid 



curve). The left panels show the results for haloes at 2: = 
while results for the z = 2 haloes are presented in the right 
panels. Top panels correspond to spin parameters computed 
using all dark matter particles within r2oo , while the bottom 
panels are for the inner region (0.25 r2oo). Table [3] lists val- 
ues for the mean, standard deviation, median and skewness 
of A for all runs, and at 2: = 0, 1 and 2. 
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Table 3. The mean, median, standard deviation and skewness of the distribution of ln(A) parameters for all haloes in each simulation 
run at z = 0, 1 and 2. We compare the distribution of spin parameters computed using dark matter particles within r200 and that 
computed using only the dark matter particles within the central 0.25r200 region. Haloes are required to contain at least 1000 particles 
within r2oo- Here In(Ao) denotes the mean, ln(Ajnt,d) the median and a the standard deviation of the distribution of InA, see equation 
Errors represent the la confidence intervals and have been determined using bootstrap resampling (1000 bootstrap realisations have 
been used). 



Simulation 



In(Ao) 



Within r2oo 
In(Amed) 
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At z = we find In(Ao) = -3.34 ± 0.03 (or Aq = 0.036) 
and cr = 0.62 ± 0.01, for our DMONLY WM^P3 run, when 
all dark matter particles within r2oo are considered. These 
values are in goo d agreement with th ose found in previous 
analyses such as iBullock et al] who obtained best- 

fit val ues of Ao = 0.035 and a = 0.5; iBailin fc Steinmet j 
ll2005l). who measu red Aq = 0.035 and a = 0.58; and 
iMaccio et all (|2008l '). who found a mean value of Aq = 0.034 
and a = 0.59. Our median spin parameter (Ame d = 0.037) 
is also in good agreement with the analysis of iBett et al.l 
l|2007l) , who found A^cd = 0.0367 - 0.0429 (depending on 
the definition of the halo). The spin parameters of the higher 
redshift {z = 1 and 2) samples are found to be slightly higher 
than those at z = 0. For example, at 2: = 2, the mean spin 
parameter of the DMONLY H^M^P3 haloes is Ao = 0.038, 
around 8 per cent higher. At least part of this (statistically 
significant) shift can be explained by a very weak depen- 
dence of spin with halo mass (see below), given that the 
mean mass of our haloes at 2; = 2 is lower than ai z = 0. 

We also checked whether excluding unrelaxed haloes 
from our sam ple made any difference to our results, as 
iMaccio et al.l (2006;) found that this reduced the mean spin 
parameter by ~ 15 per cent and lJeeson-Daniel et al.l (|201l|) 
found that spin correlates with the relaxedness of the halqj 
(they demonstrated that this is due to the strong anti- 
correlation of both spin and relaxedness with concentration). 



^ A low value of relaxedness corresponds to a relaxed halo. 



Our results are broadly in agreement: at 2 = we find a re- 
duction in Ao of around 12 per cent while at z = 2, this 
increases to 17 per cent. 

At z — 0, we find that the inner region (r < 0.25 r2oo, 
shown in the bottom left panel of Fig. [1} is characterised by 
nearly the same distribution as the overall halo. However, 
this does not hold at 2: = 2 (bottom right panel) where we 
find that the mean spin parameter is significantly higher 
(In(Ao) = -3.10 ± 0.01) than that of the whole halo (In(Ao) 
= — 3.26±0.01). Interestingly, this difference is even stronger 
for the relaxed sample at 2: = 2, where the inner region of the 
relaxed dark matter haloes exhibits a mean spin, In(Ao) = 
—3.21 ± 0.01 while the mean spin computed over the whole 
halo is found to be In(Ao) = -3.42 ± 0.01 . 

We have also considered the effects of cosmology on the 
spin parameter distribution of dark matter haloes, compar- 
ing the spin distributions from dark matter only simulations 
run with the WMAPl, WMAP3 and WMAP5 cosmological 
parameters at z = 0. The mean, median, standard deviation 
and skewness for the distributions of the log of the spin pa- 
rameters ar e also listed in Tab le |31 It is evident from this 
table (as in IMaccio et al.ll2008l ) that the spin distribution is 
not sensitive to the exact choice of cosmological parameters. 
We note that the largest change in the models we have con- 
sidered is the value of erg, which varies from 0.7 in WMAP3 
to 0.9 in WMAPl. Haloes in a lower erg cosmology are ex- 
pected to form later an d be less concentrate d than those in 
a higher erg cosmology (|Navarro et al.|[l997l 'l. The spin dis- 
tribution itself remains almost unchanged for the range of 
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Figure 2. The left (right) panel shows the spin parameter versus M200 in the dark matter only simulations at z = (2). In these panels 
we show the median value within each mass bin and errors represent the Icr scatter in the mass bin. Haloes are required to contain at 
least 1000 particles and mass bins at least 20 haloes. We plot the results from the three cosmologies we have considered: WMAPl is 
depicted using red squares, WMAP3 using green circles and WMAP5 using blue triangles. For clarity we show the median of the whole 
WMAP3 sample as a solid line and present this value in each panel. We note that the spin parameter remains unchanged as we vary the 
cosmological parameters and that there is also no clear dependence on mass. 
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Figure 3. This figure illustrates how the DM halo spin depends on the central baryon fraction of the halo (used as a proxy for galaxy 
formation efficiency). The z = (2) haloes are shown in the left (right) panel. We plot the median values and error bars show the Icr 
halo-to-halo scatter. The DM spin parameters are computed within 0.25 r2oo while the baryon fraction is calculated within 0.05 r2oo- 
The two quantities are positively correlated at z = but the correlation is absent at 2; = 2, where only the run without feedback 
(NOSN_NOZCOOL) shows a significant increase in the median A. 



parameters explored here, although we note a slight decrease 
in the spin parameter for the WMAPl haloes, consistent 
with the anti-correlation of s pin with concentration (as dis - 
cussed in lMacci"F et al.ll2008l and|j eeson- Daniel et allbOllf ). 
Haloes formed in this cosmology are expected to have the 
highest concentrations and hence may be expected to have 
the lowest spin values. 

We explore the relationship between spin and mass in 
Fig. [21 In this figure we show the median spin value and 
the la halo-to-halo scatter for several mass bins, for each 
of the three cosmologies that we have considered at z = 
(left panel) and for the WMAP3 cosmology only a,t z — 2 
(right panel) . We use least-squares fitting of the median in 
each mass bin to determine the slope and intercept of the 
A — logjQ(M2oo) relation in each simulation set, assuming 
a pivot mass of Afpivot = 10^^ M©. It is clear that the 



spin parameter is insensitive to halo mass; a weak trend is 
seen (which may explain the slight redshift dependence seen 
earlier, due to our two samples covering a different mass 
range), but is not statistically significant (we find the best- 
fitting slope to be consistent with zero at the Icr level in all 
cases - where the Icr confidence is determined using boot- 
strap resampling). Once again, we can clearly see that the 
effect of cosmology on the spin parameter is negligible. 



3.1.2 Effect of baryons on DM halo spin 

Following on from the dark matter only case, we now con- 
sider the distribution of DM halo spin parameters computed 
from the baryon runs. We only include dark matter par- 
ticles in our determination of A to directly assess the ef- 
fects of baryons on the dark matter component. The mean, 
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median, standard deviation and skewness for the distribu- 
tion of spin parameters in each run is presented in Table [S] 
For ease of comparison with the dark matter only case, we 
show in Fig. [T] the mean value for two of the most extreme 
feedback schemes we have considered: the no feedback case 
(NOSNJ^OZCOOL; green arrow) and the stellar and AGN 
feedback run (AGN; red arrow). 



From Fig. [T] and Table [3] we see that, at 2: = 0, the 
baryons induce a small increase in the mean (and median) 
value of A within r2oo. For NOSNJSfOZCOOL (where feed- 
back is absent), the increase is around 15 per cent, but for 
the more realisti(13 AGN run (where feedback is strongest), 
it is only 5 per cent. At higher redshift {z = 1 and 2), the 
effect of baryons on the mean spin parameter is weak or 
absent in all cases except for NOSNJN'OZCOOL. It is also 
interesting to note that, within the central region (0.25r2oo)) 
the baryon runs have distributions that are skewed to lower 
values than the DMONLY run (also true for the whole halo 
at z = 0) . It is possible that this is a result of the transfer of 
a ngular momentum fr om infalling satellites. As in the model 
of lMaller etlo] l|2002h . if gas is able to cool and form a dense 
clump, it will lose angular momentum via dynamical friction 
and will play an important role in shaping the lower end of 
the angular momentum distribution. 



The effects of baryons are stronger in runs with weak or 
no feedback, when we consider dark matter parti cles within 
the c entral region (consistent with the findings of lBett et al.l 
bold , who studied the effects of baryons using a more re- 
stricted set of simulations). There, the mean spin parameter 
increases by 50 per cent in NOSN3JOZCOOL (but again, 
is almost unchanged in AGN). This effect is shown more 
clearly in Fig. [31 where we plollj the median spin parameter 
for the central region (0.25 r2oo) against the baryon frac- 
tion within the very centre (0.05 rgoo)- The latter is used as 
a proxy for galaxy formation efficiency, as weaker feedback 
leads to stronger gas cooling and a larger central baryon frac- 
tion, associated with the galaxy at the centre of the halo. It 
is clear that there is a positive correlation between the two 
quantities at 2 = but this correlation has all but disap- 
peared at 2 = 2. 
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Figure 4. The ratio of specific angular momentum of the dark 
matter within 0.25r200 from the baryon runs j'dmg (< 0.25r20o) 
to that of the corresponding dark matter only run within a ra- 
dius containing an equivalent mass (jdmo(< ^0)). This ratio is 
shown against the DM mass ratio within 0.25r2oo for the 500 
most massive haloes that have been matched between runs at 
^ = 0. Symbols depict the median values while the error bars 
represent the Icr scatter. 



3.1.3 Origin of excess spin in the DM 

From equation[T] it is clear that the increase in A must be due 
to an increase of the specific angular momentum of the dark 
matter, or a decr ease in the enclosed mass. As discussed in 
iBett et al.l(|2010l '). an increase in specific angular momentum 
can result from two factors. The first is the transfer of an- 
gular momentum from gas to dark matter via tidal torques 
and dynamical friction. The second is the contraction of the 
dark matter in response to the deepened potential well of 
the system (assuming angular momentum is conserved). 

To ex p lore this further, we apply the approach of 
iBett et all (|2010l ) to our haloes at 2 = 0. We begin by 
matching haloes between the DMONLY and baryon runs; 
a match is identified by selecting the halo from the baryon 
run that contains the most of the 1000 most-bound parti- 
cles from the dark matter simulation. We then compute the 
DM mass ratio (within 0.25 r2oo) for each matched pair of 
haloes. If the DM halo has contracted due to the presence of 
baryons, the mass should increase with respect to the mass 
in the DMONLY halo. We then compare the specific an- 
gular momentum of the dark matter within 0.25 r2oo from 
the baryon run jdmg(< 0.25 r2oo), to the corresponding halo 
in DMONLY, within a radius ro containing the equivalent 
mass, jdmo(< ro). If the halo contracts due to the presence 



We find that baryons do not influence the mass depen- 
dence of the spin parameter (in almost all cases, the slope 
is consistent with zero at la and in all cases it is consistent 
with zero at 2a). The trends found in this analysis are also 
preserved when only the sub-set of relaxed haloes is consid- 
ered. 



^ The AGN run reproduces the z = observed relations be- 
tween black holes and the mass and velocity dispersion of their 
host galaxies l lBooth fc Schavc 2009). as well as the observed op- 
tical and X-ray prope rties of the groups in which they reside 
llMcCarthv et al.|[201ol) . 

^ Note: To ensure that mass dependent trends do not influence 
our comparison we do not include the 400/i~^ Mpc box dark mat- 
ter simulation when directly comparing runs at z = 0. This leaves 
us with a sample of 4329 haloes in the DMONLY simulation at 
2 = with a median mass of 8.3 X lO^'^ Mq comparable to 
the median mass in the baryon runs. 
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of baryons and conserves its angular momentum (i.e. if there 
is no transfer from the baryons) then this ratio will be unity. 

The results for each baryon run are shown in Fig. IH 
where we show the ratio of specific angular momentum ver- 
sus the mass ratio within 0.25 r2oo for the 500 most massive 
haloes that we are able to matclfl between runs at z = 0. 
From this figure, it can clearly be seen that the AGN feed- 
back acts to reduce the mass in the central regions, but as 
galaxy formation efficiency increases, we see the effect of 
halo contraction. T his effect was stud ied in detail for the 
same set of runs bv lDuffvet ail (|2010l ). 

If we now consider the ratio of specific angular momen- 
tum, we find that the AGN result is consistent with no trans- 
fer, but as galaxy formation (and thus halo contraction) is 
more efficient, a net transfer from t he baryons to the dark 
matter does occur (this is also seen in lKaufmann et"al ]|2007l : 
lAbadi et al1l2010l : ISett et al.ll201Gl '). In summary, the small 
increase in spin parameter (within 0.25 r2oo) in the AGN run 
appears to be driven by a slight decrease in halo mass, but 
in all other runs with more efficient galaxy formation, the 
increase is due to a net transfer of angular momentum from 
the baryons. 

3.2 Halo Shapes 

3.2.1 Dark matter only simulations 

Before turning our attention to the effect of baryons on halo 
shape, we will analyse the DMONLY runs. Of particular 
interest is the relationship between the shape parameters, 
e = b/a, s = c/a and T — {a? — h^)/{a^ — (?), and mass, 
M200, of the dark matter halo. Previous work has shown 
that e and s are negati vely correlated w i th mass, and T pos- 
itively correlated (e.g. iBett et aLll2007l : iMaccio et al]l2008l : 
Ijeeson-Daniel et al.ll201lh . Our dark matter only results ex- 
tend this work by studying the dependence of these correla- 
tions on both cosmology and redshift, as well as establishing 
the baseline for the runs with baryons. 

The mass dependence of the shape parameters is shown 
in Fig.O In each panel we show the median value for e, s and 
T, within several mass bins, while the error bars represent 
the la halo-to-halo scatter. Again, haloes are required to 
contain at least 1000 particles and the mass bins at least 20 
haloes. We use least-squares fitting to determine the slope 
and intercept of the mass relation for each parameter, us- 
ing a pivot mass of Mpivot ~ 10^^ M©. In the left pan- 
els we show the median values for all three cosmological 
models at 2; = 0, but for clarity, only plot the fit to the 
DMONLY WMAP3 haloes. The corresponding results for 
the DMONLY WMAP3 run at 2 = 2 are shown in the right 
panels. The best-fit parameters for the mass dependence of 
all simulation sets are presented in Table [3] 

As with the spin parameter, it is clear from the left 
panels of Fig. [S] that there is no significant difference be- 
tween the different cosmological models studied here. At 
2 = there is a strong trend for more massive haloes to 
have smaller axis ratios. Quantitatively, the sphericity varies 
from ~ 0.6 for haloes of mass 10^^ h-^ Mq, to ~ 0.5 for 

A match is found in 96% of the cases (and at least 85% of these 
matches have 771200 ratios within 20% and minimum potential 
positions within 20% of r200 of each other). 



haloes of mass > lO'^* Mq. Since e = b/a decreases 
faster with mass than s = c/a does, the triaxiality parameter 
increases with mass; over the same range in mass, the triaxi- 
ality parameter increases from ~ 0.6 to ~ 0.8. In short, high- 
mass haloes are less spherical and more prolate than their 
lower-mass counterparts. This is likely a conse quence of their 
more recent formation time (as discussed in ISpringel et al.l 
'2004 'Jeeson-Danie l et aLllioiH). Our results are in agree- 
ment with previous studies |Bett et al...200^ : iMaccio et al.l 
I2OO8I : IJeeson-Daniel et aLlbOllI ). 

At 2 = 2, we find that the correlations between the 
shape parameters and mass are much weaker, but still 
present. For example, s oc M^qo"^^ at 2 = but s oc M2oo'°^^ 
at 2 = 2. A similar trend i s seen in the evo l ution of the 
conce ntration-mass relation (|Zhao et al. 1 12OO3I : iDuffv et al.l 
I2OO8I ) and again, is likely to reflect that haloes have had less 
time to collapse, so exhibit a narrower range in formation 
times (and therefore internal p roperties). This explanatio n 
is consistent with the findings of lJeeson-Daniel et al.l l|201lh , 
who used a principal components analysis to show that con- 
centration, which they found to be essentially equivalent to 
formation time, is the most fundamental property of dark 
matter haloes in that scatter in this quantity accounts for 
much of the scatter in the other dimensionless properties of 
dark matter haloes, including the shape parameters investi- 
gated here. 

3.2.2 Effect of baryons on DM halo shape 

Baryons are known to have a significant effect on the shapes 
of haloes. As gas cools and condenses within the centre of 
a halo, it reduces the fraction of box orbits present in the 
dark matter and stars, resulting in a more spherical struc- 
ture. Feedback prevents cooling and acts to expel gas from 
the central region, resulting in a lower central mass concen- 
tration, more box orbits and a more triaxial shape (see, for 
example, iDebattista et all 1200^ : iBrvan et aLll2012f ). Here, 
we quantify the impact of these effects for our range of feed- 
back prescriptions, over a wide range in halo mass, and at 
low and high redshift. 

The mass dependence of the halo shape parameters in 
the baryon runs is summarised in Table |31 Again, this is 
for the shape determined within 7-200 using only the dark 
matter particles. From this we see that the trends seen in 
the DMONLY runs are still present in the runs with baryons. 
However, it is interesting to note that the baryons are able to 
affect the shape of the halo out to r2oo, albeit in a minor way. 
The axis ratio intercepts (both b/a and c/a) systematically 
increase with increasing galaxy formation efficiency, while 
the triaxiality parameter decreases accordingly. 

The baryons have a larger impact on the shape of the 
central region of the dark matter halo. To quantify this, we 
compare the median values for e, s and T within 0.25r2oo 
to the central baryon fraction (within 0.05r2oo) for each of 
the simulation runs. Note: To ensure that mass dependent 
trends do not influence our comparison we do not include 
the 400h~^ Mpc box dark matter simulation at 2 = when 
directly comparing runs. These results are presented in Fig. 
[6] From this flgure we can see that there is a clear trend for e 
and s to increase with increasing baryon fraction, while the 
triaxiality decreases with increasing baryon fraction. This 
trend is apparent both at 2 = and 2 = 2. At 2 = 0, 
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Figure 5. The panels show, from top to bottom, how the axis ratios (e = b/a, s = c/a) and triaxiality (T) of dark matter haloes in the 
DMONLY runs scale with halo mass (M200) at 2 = (left) and at z = 2 (right). We show the median value within each mass bin and 
the error bars represent the la intrinsic scatter. The results from the WMAPl haloes are depicted using red squares, WMAP3 using 
green circles and WMAP5 using blue triangles. Haloes are required to contain at least 1000 particles and mass bins at least 20 haloes. 
Least-squares lines of best-fit (assuming a pivot mass of 10^^ a-re computed for each simulation and are presented in Table|4l 

The best-fit line for the WMAP3 haloes is shown in each panel. 



the median sphericity is increased by ~ 25 per cent when 
comparing the extreme (DMONLY and NOSN JslOZCOOL) 
cases, while the triaxiality is reduced by approximately the 
same factor. The change is still significant in the more real- 
istic AGN model, where the sphericity increases by around 
10-15 per cent and the triaxiality decreases by a similar fac- 
tor. As expected, these results are consistent with the idea 
that gas cooling to the centre of haloes results in more spher- 



ical haloes. It also drives down the value of the triaxiality 
parameter, tending to make the haloes considerably less pro- 
late than seen in dark matter only runs. Feedback reduces 
the amount of gas that is able to cool and condense to the 
centre of the halo (resulting in a lower central baryon frac- 
tion) and hence reduces the impact on the shape of the halo. 
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Table 4. Best- fit slope and intercept values for the mass dependence of halo shape parameters (e, s and T), assuming a linear relationship 
between each shape parameter and logjQ(M20o)i about a pivot mass of lO^'^ /i^^Mq. Fits are performed for dark matter particles within 
7'200i for all runs considered at z = 0, 1 and 2. The errors correspond to Icr confidence intervals, determined by bootstrap resampling 
1000 realisations of the data. 
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4 SUMMARY 

In this paper, we have exploited a subset of r uns from the 
Overwhelmingly Large Simulations (OWLS; ISchave et al.l 
l201Cf ) to investigate the impact of baryons (through gas cool- 
ing, star formation and feedback from stars and black holes) 
on the spin and shape of dark matter haloes. Our results al- 
low statistically meaningful conclusions to be drawn regard- 
ing the impact of baryons on these properties, due to the 
large number of haloes spanning a wide dynamic range in 
mass. We have also checked whether our results depend on 
cosmology and redshift. Our main results are summarised 
below. 

• The spin distribution of dark matter haloes in simula- 
tions without baryons is characterised by a log-normal curve, 
with best-fit values of Ao = 0.036 (0.038) and a = 0.62 

1 0.60) at 2 = ( 2 ), in agreement with previous work 
Bullock et all I2OOII : iBailin fc Steinmet j l2005l : iBett et all 



I2OO7I : iMaccio et al.ll2008l ). The distribution is very similar 



for the WMAPl, WMAP3 and WMAP5 cosmologies, sug- 
gesting that there is no strong dependence on as (the pa- 
rameter that varies the most between the three models) . No 
significant dependence of spin with mass is seen, both at 
z = and z = 2. At z = the spin parameter remains es- 
sentially unchanged if computed using o nly mass within the 
central region (0.25r2oo), as found by iBailin fc Steinmeta 
l|2005h . However, at 2 = 2 the inner region of haloes has a 
higher mean spin than that computed over the whole halo. 
Restricting our sample to relaxed haloes causes a small (10- 
15 per cent) decrea se i n the mean value (in agree ment with 
iMaccio et al1l2006l and IJeeson-Daniel et al1l2011^ . 



• The spin distribution of dark matter haloes extracted 
from the baryon runs is not significantly different to that 
of dark matter only haloes when computed using all dark 
matter particles within r2oo- However, in the central regions 
(0.25r2oo), where baryons are expected to play an important 
role, haloes in runs with absent or weak stellar feedback tend 
to have higher median spin values than those from stronger 
feedback runs (which are very similar to the dark matter 
only case). We showed that this is, at least in part, due to 
the transfer of angular momentum from the baryons to the 
dark matter in the former runs. 

• Dark matter only haloes extracted from OWLS typi- 
cally have sphericities of ~ 0.5 to 0.6 and triaxialities of 
between 0.6 and 0.8 (indicating triaxial to prolate shapes) 
over the mass and redshift ranges we have explored. More 
massive haloes have more spherical and less prolate shapes. 
Again, we find that halo shape is insensitive to the choice of 
cosmological model. Galaxies and groups at z — 2 show the 
same trends with mass as the groups and clusters at z — 0, 
but weaker. 

• When baryons are included, we find that the mass de- 
pendent trends remain, and that the intercepts of the re- 
lation between sphericity (triaxiality) and mass slowly in- 
crease (decrease) with increasing galaxy formation efficiency. 
At M200 = lO^^/i~^M0, baryons increase the dark matter 
shape parameters by around 10 per cent in the most ex- 
treme case (no feedback). A similar result is seen at higher 
redshift. Larger differences are again seen when we consider 
only the central regions of the halo. 
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Figure 6. The dependence of the shape of the central halo on galaxy formation efficiency (as measured by the baryon fraction in the 
inner region of the halo). Here we show how the halo shape (computed within 0.25r2oo) depends on the baryon fraction contained within 
0.05r2oo- We show the median and Icr halo-to-halo scatter for the each of the feedback implementations. All haloes contain at least 1000 
particles. Results at z = (2) are shown in the left (right) panels. Runs with more efficient galaxy formation (weaker feedback) contain 
haloes that are more spherical and less prolate than runs with stronger feedback. 



In conclusion, we find that the baryons have a very mi- 
nor efTect on the spin and overall shape of the entire dark 
matter halo when the feedback is strong enough to match 
observed stellar mass fractions. In particular, the model with 
AGN feedback can reproduce several observational proper- 
ties on galactic and groups scale at z = by removing gas, 
suppressing the baryonic impact on the dark matter halo 
shape. It should therefore be safe to assume results from 
dark matter only simulations when considering the overall 



halo properties, at least on the scales resolved by our simu- 
lations. However, even when feedback is strong we find that 
baryons have a significant effect on the shape of the inner 
halo. 
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APPENDIX A: RESOLUTION TESTS 

In order to quantify the effects of resolution, the shapes of 
haloes extracted from a 512'' particle run with a maximum 
force softening length of 2 ft" ^ kpc (used in this analysis) are 
compared with those from a corresponding lower-resolution 
run (containing 256^ particles with a maximum softening 
length of 4 kpc). Resolution tests for the central re- 
gions (0.25r2oo) of dark matter only haloes and the haloes 
from the weak stellar feedback run are shown in Fig. \M 
while the resolution tests for the overall properties com- 
puted within r2oo are shown in Fig. IA2I These figures show 
the spin (A), e = b/a, s = c/a and triaxiality (T), respec- 
tively, as a function of halo mass. In each of these plots the 
properties of haloes extracted from the 512"' particle simu- 
lations are shown with filled symbols and the results from 
the lower-resolution run with 256^ particles are shown with 
open symbols. Dark matter only runs are shown as blue 
squares while the weak stellar feedback runs are shown as 
black circles. The error bars are the one-sigma bootstrap re- 
sampled median distributions (1000 bootstrap samples have 
been used). Vertical lines show the 1000-particle cuts that 
have been used in this analysis. Clearly, the properties we 
discuss are well resolved beyond 1000 particles. 
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Figure Al. Resolution tests for halo properties determined within the central 0.25r2oo- Properties of haloes extracted from the 512"^ ( 
256^) particle simulations are shown as filled (open) symbols. Dark matter only (weak stellar feedback) runs are shown as blue squares 
(black circles). The error bars show the Icr bootstrap resampled median distributions (1000 samples). We note that above 1000 particles 
(vertical lines), all quantities are well converged with respect to the numerical resolution, and we therefore adopt this resolution limit. 
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Figure A2. As above but for halo properties computed within r2oo- 



